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The Gross-Pitaevskii equation is a widely used model in physics, in particular 
in the context of Bose-Einstein condensates. However, it only takes into account 
local interactions between particles. This paper demonstrates the validity of using 
a nonlocal formulation as a generalization of the local model. In particular, the 
paper demonstrates that the solution of the nonlocal model approaches in norm the 
solution of the local model as the nonlocal model approaches the local model. The 
nonlocality and potential used for the Gross-Pitaevskii equation are quite general, 
thus this paper shows that one can easily add nonlocal effects to interesting classes 
of Bose-Einstein condensate models. Based on a particular choice of potential for 
the nonlocal Gross-Pitaevskii equation, we establish the orbital stability of a class 
of parameter-dependent solutions to the nonlocal problem for certain parameter 
regimes. Numerical results corroborate the analytical stability results and lead to 
predictions about the stability of the class of solutions for parameter values outside 
of the purview of the theory established in this paper. 
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I. 



INTRODUCTION 



The last 15 years has a seen a rapid growth in interest concerning the modeling of Bose- 
Einstein condensates. The body of literature concerning this subject is too vast to consider 
here, but a simplified description of the field would include the study of the Gross-Pitaevskii 
equation 



where a — ±1, with +1 corresponding to repulsive interactions between particles in the 
condensate, and —1 corresponding to attractive interactions. The function ip represents an 
approximation to the wave function used to describe the probability density for the location 
of particles in the condensate. 

The validity of this equation as an approximation to the many-particle formulation of 
the problem has been established in pp. However, an assumption of a pairwise 5-function 
interaction among particles is used to derive Q. This clearly cannot capture all of the 
physics in the problem since each particle in the condensate exerts forces that act at a 
distance. Thus the next order of approximation to the many-particle formulation would 
be to include a more general interaction potential simulating nonlocal interactions between 
particles. This is done in j2] by studying the modified one-dimensional Gross-Pitaevskii 
equation 



where e > 0, and R(x; e) = ~C(~), with ((x) being a positive, even function such that 



in the sense of distributions. In [2], e is called the nonlocality parameter. The authors of 
[2] assume that the condensate is trapped in both a harmonic confining potential and an 
external standing-wave potential. While in [2] a three-dimensional version of ^ is derived, 
the presence of the standing-wave potential allows the reduction to a one-dimensional model 



Nonlocal models like (J2J) are also called Hartree-Fock equations. These have been ex- 
tensively studied in the case that ((x) = l/\x\, i.e. in modeling Coulombic interactions 
between particles (cf. jl], [5]). Recent literature on the formation of dipolar condensates 
has introduced nonsingular nonlocalities characterized by cubic decay (cf. j6], [7]). These 




(1) 





lim R(x; e) = 5(x) 



(cf. i). 
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nonlocalities with cubic decay fit into the class studied in this paper. Other models with 
a varying nonlocality parameter have appeared in the optics literature [8]. The analysis of 
the well-posedness and convergence of nonlocal, nonlinear Schrodinger type models to local 
ones can be found in [9] and |10| . though the models examined in those papers are different 
from those studied in this paper. 

The authors of working with the potential V(x) = Vq sin 2 (kx), derived the traveling- 
wave solutions 

^{x,t)=r sol {x)e w ^\ (3) 



where 



rL(x)= B--^sin a (fcr), 
tan( W )=Jl-^|_tan(M 



V + k 2 : „ Vq 



u (k) — + aB 



o 



2/3(k;e) 

f3(k;e) = I R(x; e) cos(2kx)dx, 



with B a constant called the offset size. Defining A = a ^.^ , the traveling- wave solution 
can be rewritten as 

ip(x,t) = (^/B cos(kx) + iV B + A sin(A;x) j e 



iuit 



which shows the spatial component of ^ is periodic with period 2n/k. The coefficients 
appearing in ip(x,t) must satisfy the restrictions 

B > max{-A 0}, a = ±1, (3(k; e) ^ 0. 

Setting a and k equal to one, and taking 5 = 1, which in [2] is described as large, and 
Vq = —1, the authors of [2] study the stability of ^ by numerical simulations using ((x) = 
e~ x2 and V(x) = Vosin 2 (fcx). The authors report results which numerically demonstrate 
that for the local case, i.e. when e = 0, (|3]) is stable with respect to perturbations due to 
roundoff error in the numerical simulation. However, their results also suggest that ^ is 
unstable when the nonlocality parameter e is positive, and that the instability emerges at a 
fixed time in their simulations, independent of the value of e. The authors of [2] also study 
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the effect of changing the convolution kernel, and they report that the results are similar to 
those for the case ((x) = e~ x2 . 

It is conjectured in [2] that a beyond- all-orders phenomena may be responsible for the 
behavior exhibited in their numerics. As pointed out in [2], if the behavior exhibited in 
their numerics is accurate and truly independent of the choice of interaction potential, then 
(|2j) cannot be viewed as a valid generalization of 0. That is to say, no matter how small 
one makes the nonlocal interaction term, the results of J2] seem to imply that one cannot 
approach the local behavior. This is described as a lack of asymptotic equivalence of stability 
(AES) in 0. 

The purpose of this paper is to address both the issue of whether or not ^ is AES to 
([!]) and under what conditions ^ is a stable solution of Equation (|2]). To do this, we first 
fix some notation and introduce the spaces in which we work. Let St denote the circle 
of circumference T. Introduce the space L 2 (St) which is the completion of the continuous 
T-periodic functions in the norm 

/ , \ 1/2 / ,T/2 \ 1 I 2 

\\f\y { s T) = \ K J s \m\ 2 dxj = U T/2 \f(*)\ 2 dx 

In practice the integral f s \f(x)\ 2 dx could be evaluated over any interval of width T since / is 
a T-periodic function. Note, throughout the remainder of the text 1 1 " 1 1 x 2 (St) * s a bbreviated by 
1 1-| | 2 . We define the norm, denoted as | j • j | 2 w , of the product space L\(St) = L 2 (St) X L 2 (St), 
via 

' ' f 
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2 , II ||2 
2 + 2 



2,v 

For operators B that map L 2 (St) to itself, we denote the norm of B via 

||£>|| 2 = sup ||jB/|| 2 . 
Il/ll 2 =i 

The norms ||jB|| 2u are defined in an identical way. The Sobolev spaces H s (St) are defined 
as follows: 

{00 
feL 2 (s T )-. ori/ii 2 <°° 
j=~oo 

where 

4n 2 j 2 



(j) = 1 



J-2 



1 



5 



and the terms f n come from the Fourier series of f(x), which is 



oo 



/(*) 



where 



1 



T 



3 = 




(4) 



1 



e -2mjx/T ■ ^ 



v 




J-T/2 



where e* denotes the complex conjugate of ej. The product space H s (St) X H s (St) is 



denoted by H^(St)- Finally, define the Fourier transform of h(x) G £ 2 (M), say h(s), by 



h(s) = / e~ iSx h(x)dx. 



J —oo 



To address the issue of AES, we first prove the local-well posedness, in H s (St) for s > 1/2, 
and global- well posedness of ^ over the space Hi(St) based on the following assumptions. 

HI: The potential V(x) is a smooth, T-periodic function, 

H2: CO) > 0, 

H3: C{x) G L\R) with \\(\\ Ll(M) = f R C(x)dx = 1, 

HA: x((x) G L\R), and 

H5: \((s)\ < (1 + l^l)- 1 / 2 ^, where e > 0. 

Note that the maximum of ( could be chosen larger than one without affecting our results. 
We make this choice in order for a cleaner presentation. Using the local and global-well 
posedness results, we prove 

Theorem 1. Let a = 1. Assuming the hypotheses HI — H5, choose constant C > such 
that Wipo^W Hl r s T ) — ^ f or e — w here ipo, e (x) is an initial condition for §2§ and ipo,o(x) is 
an initial condition for Q. Let ^o,e — > ^0,0 in the H\(St) norm as e — > + . Let ip(x,t) 
and ip e (x,t) be the unique T-periodic solutions to ([!]) and ^ respectively for the given 
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initial conditions. Then there exists a constant C u > and a function g(T;e,e), where 
lim^e' q(T; e, e) — for e, e' > ; such that, for any finite T > 0, we have the bound 

\\M;t)-M;t)\\L«>(s T ) < (\\^o,e-^o,o\\ HliST) + SC u fg(f;e,0)^ e 3C ' f 
forte[0,f\. 

Thus, on any finite interval of time, as one lets the nonlocality parameter approach zero, 
the solution to ^ converges uniformly in space to the solution of ([lp. This shows that 
AES is a common feature for a large class of potentials and nonlocal, repulsive interactions. 
Therefore the results in [2] are likely due to artifacts of their numerical computations, as 
opposed to being inherent to the equation. 

As to the stability of (|3]), we first need to define the notion of stability to be established 
(cf. [llj). Let ?p(x,t) denote a solution to either Q or Q with initial condition ip(x, 0). 
Writing ^ as (p UJ (x)e~' iU!t , we say that ^ is orbitally stable in Hi(St), if for any p > 0, 
there is a 5 > such that if H^oO^) — ^{ x )\\hi(S t ) < ^ ^ en 

sup inf - (j)Jx)e lc \\„ ,„ > < p. 

The other notion of stability we use is that of spectral stability. First, separate Q or ^ 
into real and imaginary parts. Denote the linearization of either of these systems around (J3]) 
as JL. Using the scaling x — > kx, JL has terms that are 27r-periodic functions. Let a(JL) 
denote the spectrum of JL computed over the space ii/f (SWn)) n £ N. In effect, we are 
computing the impact of perturbing ^ by 27rn-periodic perturbations, or 27rn/A;-periodic 
perturbations in the unsealed coordinate. We say ^ is spectrally stable if for A G a(JL), 
Re(A) < 0. Note, more details are provided in Section 3. Also, given that the nonlinear 
problem is Hamiltonian, the condition of spectral stability reduces to having spectrum only 
on the imaginary line, i.e. Re(A) = 0. With these definitions in hand, we prove the following 
three theorems. Throughout these remaining theorems we assume that 

Hi': V(x) = 1/ sin 2 (fca;), 

H2': C(x) > 0, C(x) is even, and ||CHli(r) = h 

H3': C > 0, and 

HA': C(s) < (1 + \S\)- 1/2 '", with e > 0. 
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Theorem 2. Let a — I. Assuming the Hypotheses Hi' — HA , for any values of k and the 
nonlocality parameter e, and for perturbations of period , where n E N, the solution ^ 
is spectrally stable for sufficiently large offset size B, V < 0, and \ V \ sufficiently small. 

Theorem 3. Let a — 1. Fix the nonlocality parameter e and the value k. Assuming the 
Hypotheses HI — HA' , for offset parameter B sufficiently large, Vq < 0, and \Vq\ sufficiently 
small, the solution ^ of the nonlocal Gross- Pitaevskii equation (|2| is orbitally stable with 
respect to perturbations with periods T = where n E N. 

Theorem 4. Let a — 1. Assuming the Hypotheses Hi' — HA' , with A = a ~^° e - ) , if A > 
2A6k 2 , then for offset size B and e sufficiently small, ^ is spectrally unstable with respect 
to perturbations of period ^r, where n E N. 

The content of these three theorems shows that the role of a small nonlocality parameter 
is dependent upon the other parameters in the problem, particularly the offset size B. 
Theorems [2] and [3] are proven by showing that if B is sufficiently large, then the operator 
L is positive semi-definite. Then, using Krein signature arguments found in [19], we get 
both spectral and orbital stability. Thus, introducing small nonlocality should not effect 
the stability of pj), while Theorem [4] shows that if B is too small, then even removing 
the nonlocality parameter does not stabilize the solution. In contrast, as is shown later, 
we can always get a spectrally stable problem by letting e — > oo for any choice of the 
other parameters. Thus it appears that while a small amount of nonlocality does not affect 
stability, large amounts do. 

The above theorems do not allow for arbitrary choices of parameters since each theorem 
requires |Vo| to be small, which ensures that L remains positive semi-definite. We cannot 
at this time provide explicit bounds on how large |Vq| can be such that Theorems [2] and [3] 
remain true since we cannot control the spectra of JL for Vq ^ 0. Therefore, we must treat 
the parameter values used in [2] as outside the scope of what is proved in this paper. We 
provide numerical experiments in order to make conjectures about the stability of ([3]). First, 
we use the above theorems to calibrate our numerics by picking parameter values that can 
reasonably be believed to satisfy the constraints of Theorems [3] and [4j Our numerics behave 
as the theory predicts. Second, we present numerical experiments using the parameter values 
found in [2], and from this we conjecture that in fact ^ should be stable for the parameter 
values chosen. As mentioned above, these are B — 1, Vq — —1, k — 1, and a — 1. 
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As in [2], a pseudo-spectral method is used for the spatial variable, while a standard 
Runge-Kutta method is used for time evolution. The most likely explanation for the discrep- 
ancy between the results reported here and those of [2] is the way in which the convolution 
is handled. In this paper, no approximation is made in the integral or to the convolution 
kernel. However, in [2] , it appears an approximation is made to the kernel which introduces 
an error that appears difficult for the pseudo-spectral method to resolve. In private com- 
munications, the authors of [2] have been made aware of these discrepancies. They have 
encouraged the explanation for them in this manuscript. 

The structure of the paper is as follows. In Section 2, we present the proof of the AES of 
([IJ and ([2]). In Section 3, we find the linearization around (|3]), and we establish some basic 
results about the convolution kernel that are used later. In Section 4, Theorems [2] and [3] are 
proved, while in Section 5, Theorem [4] is proved. Finally, Section 6 presents the numerical 
results. 



II. ASYMPTOTIC EQUIVALENCE OF STABILITY 

We proceed in the following fashion. In order to make the presentation self-contained, 
we first establish the local-in-time well-posedness of the nonlocal Gross-Pitaevskii equation 
from which we obtain a local-in-time form of AES. We then establish the continuity in e of 
solutions to ([2]). Finally, we establish the global-in-time well posedness of ^ which allows 
us to prove Theorem [TJ Note, we use Hypotheses HI — if 5 throughout the remainder of the 
section. 

We begin by establishing some basic lemmas concerning the convolution kernel R. We 
show, using the assumptions stated for Theorem [TJ that the Fourier transform of the con- 
volution kernel is Lipschitz continuous in e. 

Lemma 1. One has 

\R(s;e)-R(s;e')\ < |e-e'||S| \\x((x)\\ L i (R) ■ 



Proof. With 



/oo 
e~ iSx R{x; e)dx, 
-00 



we have 

/oo 
e- ieSx C(x)dx. 
-oo 

Using the Mean- Value Theorem, one gets 

\R(s;e)-R(s;e')\ < \e - e'\\s\ |K(-)|| il(R) • 

Thus the result is shown. Using Hypothesis if 4, which amounts to assuming H-COHl^r) < 
oo, we also have that the bound is meaningful. □ 

Let R e (-) = R(-,e). We then show 

Lemma 2. Given Hypotheses H3 and H5, for f G L 2 (St), 

R,*f = R.*(jr fMx)) = jr /a f 2 ^-) ej {x). 

\j=—oo / j=—oo ^ ' 

Proof. First we note that 

E i/iii^^i<n/ii 2 (E ^ 2 (t0) <oo > 

j=~oo \j=-oo ^ ' / 

since for j ^ 



R e *ej(x)= / R(y — x;e) j= — dy 



-oo v T 

/oo 
R(y; e)e-*^ T dy, 
-oo 



R e * e (x) = R e (0)e (x), and R e (s) < (1 + e|s|)~ 1/2 ~ e by Hypothesis H5. By Hypothesis 
if 3, \\Re(x)\\ L uiBL) = 1 5 an d by a corollary to the Dominated Convergence Theorem ([T2]. 
Theorem 2.25) one has 

/ oo \ oo 

R e*\ E A e i( x ) = E fj R e*ej(x), 



or 



V7=-oo / j=-oo 



iW= E ^(y) e i(4 



j=-oo 

and the result is shown. □ 



From the previous lemma, one gets 
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Lemma 3. Let s > |. One has that 

\\ R ^\f\ 2 \\ Hsi s T) <\\f\\k(s T) - 



Proof. By definition, one has 
\\R e *\f\ 2 \ 



A 2 



3=-oo 



Let |/| 2 = J2j qjej(x), where ej(x) is as in Q, so that the Fourier series of R e * \ f\ 2 , using 
Lemma [2j is 

'2n j 



j=-oo 



Hence, 



and therefore 



Re*\f\ 2 (x)= * A 
o 

2ttj 



(5) 



< 



\R**\f\ 2 \\ Ha{ST) < 



\H s (St) 



< 



2\A| 



ft (St) 



where the last inequality comes from the fact that H s (St) is an algebra for s > 1/2 [12J. □ 



A. Local-in-Time Well Posedness of the Nonlocal Gross-Pitaevskii Equation and a 

Weak AES Theorem 

For the small-time argument, with initial condition ip 0te (x), we rewrite ^ in the Duhamel 
form 

rt 

il> e (x,t) = e- jifla Vo,e(a;) - ia / e _tI "" ( *-' ty e (a;, t')R e * \ip e (x,t')\ 2 dt , (6) 

Jo 

where L sa = —\d 2 + V(x), with V(x) assumed to be, by Hypothesis HI, a smooth T- 
periodic function. As for controlling e~ lLsat , since the operator —iL sa is skew-adjoint, by 
Stone's theorem [13], e~ lLsat is a unitary operator from L 2 (St) to itself. One also has that 

where L sa (j) denotes the symbol of L sa . Since L sa is self-adjoint, L sa (j) is strictly real, and 
so one has the bound 

|| e iL„Y|| ^ < ||/|Ih s ( St ) , 
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since \e tLsa ^* | = 1 for all j. 

Throughout the remainder of the section, we assume L sa is acting on the space H s (St) 
where s > 1/2. Using Lemma [3j and defining G e (ip) by, 



one has for s > 1/2 

sup | |G e (V>) | | Hs(5t) < \\ipo,e\\ H Js T ) + T [ sup ||^||^ (5t) 

ie[0,f] \tG[0,f] 

We then choose constant C > such that 

ll^o, e || Hs (s T ) ^ 

for e > 0, with C independent of e. We define the metric space B r , for r > C, by 



B r = { G L^H^St); [0,T]) : sup ||V|Ih3 ( s t ) < r 

{ te[o,f] 

where 



L°°(F s (5 r ); [0,T]) = <^ ^(x,t) G H 8 (S T ) Vt G [0,T] : sup |MI H , (SW < oo 

[ te[o,T] 

The map G e takes B r to -B r for T < ^f-. It is straightforward to show, again using Lemma 
3j that G t is a contraction for f < and therefore, using the Banach Fixed Point Theorem 
[14"] . one has local well posedness for initial condition ||V ; o,e(a ; )||^ s (5 T ) < C, s > 1/2, on the 
space _B r for 

, r-C 1 

T < mm 



r 3 ' 3r 2 

From the local well-posedness result, we now prove 

Lemma 4. One has that \\i^ e ( 4 , t)\\HJS T )> ll^ e (''^)lli oo (S T )' an< ^ HV'eC'j ^)llz°°(sv) are continu- 
ous functions of time for t G [0, T] . Further, the Fourier coefficients of ip e (x, t) and ip^(x, t) 
are continuous functions of time for t G [0, T]. 

Proof. Let t, t G [0,T], and e > 0. Letting ip e (x,t) denote the solution for initial condition 
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ipo,e{%), we have 

|U e (-,t)-<M-,*)| 



Iff* (St) 



< 



-iL 3a (t—t 



H S {S T ) 



Hs(St) 



+ 



H S (S T ) 



where 

AT(^)=^(x,t') J R £ *|^(-,t')| 2 (x). 
For the first term, we have that 



O.c 



# S (S T ) 



oo 

< £ or 

j=-oo 



e ~iL sa (j)t _ e ~iL sa (j)t 



^0,e(j) 



Using the Dominated Convergence Theorem shows that this term then vanishes as t — > t or 
vice versa. The local well posedness result and Lemma [3] ensures that 

\wm\ Hs{ST) <r\ 



so we have 



-iL sa {t-t ) 



Af(ip e )dt 



< t - t 



H S (S T ) 



Using a dominated convergence argument shows that the remaining term must also vanish as 
t — > t. Thus we have shown that 1 1 -0 e ( • , ^ ) 1 1 ^ s ( s-r) ^ s a continuous function in t for t £ [0, T]. 
By a Sobolev embedding [12], we also have that \\ipe{-,t)\\ Loo , s \ * s continuous in t and 
bounded above by r. Given that 



L°°(S T ) 



<2r|hk(-,*)-^(-,t)| 



\L°°(S T ) 



we then see that HV'eC') 0IIl°°(5 t ) i s continuous as well. This result then immediately gives 
that the Fourier coefficients of ip e (x,t) and ip^(x,t) are also continuous in time since 







< < 




;t) 


V'eC'j*) | Il°°(S t ) , 






< < 




;t) 


~ 4>e(- >*) L°°(5 T ) ' 



for j ^ 0. The j = case is treated identically, so the result is proved. 



□ 
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Taking e and e to be nonnegative, we now choose the initial conditions to be continuous 
in e with respect to the x-norm, i.e. 

hmJ|V>o, e -VvlL ( s T ) = °- ( ? ) 

We then prove 

Lemma 5. For e > 0, t G [0,T], and s > 1/2, if the initial condition ip 0>e is continuous in 
e with respect to the \ \-\ \H s (s T )' norm > ^he solution ipe{x,t) is continuous in e with respect to 
the \ \-\\H a {s T )' norm ' 

lim||^(-,t)-^(-,t)||^ (5T) = 0. 
Proof. Choosing e, e > 0, one has that 

\\A'(;t)-M;t)\\ HsiST) < ||w -^A\h 3 (St) 



+ / \\A' R e' * l^e'f ~ ARe* \A \ 2 \ \ Hs {St) dt . 
J 



In the last term of the above inequality, the convolutions depend on the different values of 
e. Thus 

ip e >R e > * |Vv| 2 - 4>A * IV'el 2 = (Vv - ip e )R(-;e) * |Vv| 2 

+ ^(R(. ;e ')-R(-;e))*\^\ 2 

+^i?(-,e)*(|^| 2 -|^| 2 )- 
For t e [0, T], using the local-in-time well posedness result, the following inequality 

\\ip £ 'R £ > * \ifj e '\ 2 -ifj e R e * |^| 2 ||^ s(5t) < 3r 2 ||^ -ip e \\ Ha[ST) 



+r\\(R,-R e )*\^\ 2 \\ m 



(St) 



holds. To control the last term in this inequality, set (as in Lemma |3 

|VV(M)| 2 = ^9;0)ej(x), 
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from which one gets that 



(i2(.;e')-^(-;e))*|^(-,*)| 2 = E^*)( A 



From the local-in-time well posedness result 



27TJ 

T 



W Ml 



# S (S T ) 



< 



so that using Lemma [TJ one gets the pointwise estimate 

lim 1^)1 



T 



T 



0. 



where the index j is arbitrary. Using Hypothesis H5, R(-;e) is uniformly bounded in j, so 
the terms 



Rl^.e 



T J \ T 

are uniformly bounded for all j . Using the Dominated Convergence Theorem, one sees that 



lim 



H 3 (S T ) 



0. 



(8) 



Then 



| |W (-,*)- ^(-, *) I \ Hsi s T ) < MT, e, e ) + 3rM ^ (•, t ) - ^(-, t ) 



where 



with 



te[o,T] 



(7(t;e,e') = (£(•; e) - «(•; e)) * 



H S {S T ) 



The term H^oe' — ^o^llij (5 ) van i snes as e — )■ e' by assumption (see ([7])). Further, since, 



as shown in Lemma |4j the coefficients qj(t) are continuous in time, this makes the term 
C(t;e,e) continuous in time since it is a uniform sum of continuous functions. Thus, the 
supremum is attained at some time t*, and since ffify holds for any t £ [0, T], one has that 



lim g(T; e, e') = 0, 



where g(T; e, e') = sup tg [ ^ C(t; e, e'). Therefore 



limi(T,e,e') = 0. 
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We know from Lemma |4j that ||^> e '(-,t) — ifj e (-,t)\\ H ^ ST ^ is a continuous function in t for 
t G [0, T]. Using Gronwall's inequality [14], one gets 

11^ 0,0- M;t)\\ Hs{ST) <A(f,e,e')e 3r2t , 
and the result is therefore proved. □ 
Since s > 1/2, it follows that 

||</v(-,t)-<M-,t)IU (5T) < \\^(;t)-M-M Hs (s T ), 

and thus 

\\^(-,t) -A(-,t)\\ L ^ ST ) < MT,e,e')e 3r2t . 

Therefore the lemma establishes that (J2| is local in time AES to ((!]). Further, once a global- 
in-time well posedness result holds for §2§ (which amounts to establishing a uniform bound 
on r for all time), the above lemma immediately furnishes a global in time AES result. 



B. Global-in-Time Well Posedness and the AES Theorem 

As established in [2], ^ has at least two conserved quantities: the L 2 (St) norm and the 
Hamiltonian 

W) = \\ (l^l 2 + 2V(x)|Vf + a\^\ 2 R e * \^\ 2 )dx. 

We choose if)(x, •) G Hi(St)- Following the argument in [T5], with R e positive by Hypothesis 
H2, a — 1, and tp e {x,t) a solution to (|2]) on time interval t G [0,T], with initial condition 
■06,0 (x) G Hi(St), then one has that 



«/ .St 1 



<2|^(^,o)| + 2 M - 2 



Using Young's inequality, one has 

|2 D ... L/. 12 J™ / IL/. ||2 



|-0e,o| #e * |-0e,o| tfe < H^olliao^) ||^,0|| L 2 (St) < \ \lpefi\\ Hl (S T ) " 

Thus the assumptions guarantee that \H{^ et o)\ < 00. Since the L 2 (SV) norm is also con- 
served, there exists a constant M such that 



sup \\M;t)\\ Hl (s T ) <Cu- (9) 
te[o,f] 
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This bound is independent of t. If we now try to iterate our local-in-time well posedness 
argument onto a time interval [T — e, T), where < e 1 is chosen so that the intervals 
[0, T) and [T — e, T) overlap, then for the new interval we may let C u take the role of the 
value C. We must work on a ball Bf with f > C u and 

= [f-C u 1 1 

T < mm <^ ^ 3 , > ■ 

Since the inequality ^ is independent of time, one can repeat the derivation of Q on the 
time interval [0, T) and obtain the same bound. Thus one can iterate the local argument 
such that the value of f need not increase, and thus the width of the new intervals can be 
set to a fixed value. This establishes for the repulsive case a global existence of solutions 
to (J5]) in Hi(St) for e > 0. As argued above, one can immediately extend the argument in 
Lemma [5] so that one has a global AES theorem. 



III. STABILITY: THE LINEARIZATION AND ITS PROPERTIES 



Having established TheoremJTJ we turn to analyzing the stability of ([3]). Note, throughout 
the remainder of the paper we assume Hypotheses Hi' — if 4' as listed in the Introduction. 
Writing ^ as ip(x,t) = (p u (x)e~ lult , and introducing the transformation r = cot, we see that 
0aj is a stationary solution of the equation 

1 f°° 

iip T = —-d^if) + aip / R(x — y; e)\^(y, t)\ 2 dy + V(x)ip — coip- 



(10) 



With ip(x,r) = u(x,r) +iv(x,r), we rewrite (10) as 



u \ I uR e * (u 2 + v 2 
J<L \ + a 
r l I \ v I \ vR e * (u 2 + v 2 




(11) 



where 



J 



L 




and L = —^d^+V(x)—u. Note, (11) is posed over Hf (S^/k), but the global-well posedness 
result established in the last section carries over without issue. We set V(x) = V sin 2 (/ex) 
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from Hypothesis HI . Letting 

u u {x) = <f) u ,r( x ) = v ^ cos(x), v u (x) = <f> u ,i(x) = VB + A sm(x) , 

and equating m = u w + ew(x, r) and v = v w + e^(a;, r), and collecting all (9(e) terms, we get 
the linearized system 



W I J < L { W \ +a{ WRe * + + 2UujRe * ( Uu)W + VujZ ' 

Z / r I V Z J \ zR c* ( U l + V l) + 2v coRe * (UuW + V u z) 

V 

With A = — — -, we have 

ap(k; e) 

Re * iul + vl) = R € *(B + Asm 2 (kx)) = B + AR e ^ 1 ~ co < 2kx ^ 



2 



Thus, since R e is an even function by Hypothesis H2' , we write 



R * (ul + vl) = B + A ( LJiM. + p( k . e ) sin 2 (fcx) | . 



and we have, for a = 1, 



L + R € *(ul + v 2 J = - 1 -(d 2 x + k 2 ). 



Introducing the transformation x — > kx, so that the potential and ^ are now 27r-periodic 
functions, and defining 

\dl + l) 



2 

o -f(^ + i), 

we rewrite the linearized system with a = 1 as 



2 / \ Z / \ V u Ru e * (u^w + V U Z 



1 / X \ 

where Rk, e (x) = Rk{x; e) = — R ; ej. Defining 



Rk 

and letting 



Rk,e* 

R k , e * 



D= + ^ (12) 
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we can rewrite the linearized system as 



with the operator L given by 

(cos(x) \ / cos(:r) Dsin(x) 
Rk,e 

-Dsinfa;) / \ cos(x) -Dsin(x) 

Using separation of variables, i.e. w(x, r) = w(x)e Xr and v(x, r) = v (x)e Xr , formally gives us 
an eigenvalue problem. We now study the spectrum of JL over H%(S 27Tn ) C L\(St)- Note, 
the fact we are working over the space i?f (SWn) reflects the fact that we have separated the 
perturbations of the exact solution into real and imaginary parts. 











JL j 






I 











A. The Eigenvalue Problem on S^n 

We wish to solve the spectral problem 

, w(x + 2rcn) = w(x), z(x + Inn) = z(x), 

where n E N. As will be shown after this section, the operator JL on the domain D(JL) = 
Hl(S2Tm) C L\(S2Tm) has a compact resolvent operator. Therefore the spectrum, a(JL), of 
the operator JL is discrete, and solving the eigenvalue problem is sufficient to determine 
the spectrum. To find the spectrum of JL, we note that an arbitrary 27m-periodic function, 
f(x), can be decomposed as 

n-l 

f(x) = Le~^ = £ Le~^ x = £ f r (x)e*°, 

m m r=0 

where f r (x) is a 27r-periodic function, and m = rmod n. Therefore, one can apply a similar 
decomposition to w and z so that 





One can show, for real //, that 

JL{e 1 ^-) = e^JLp 
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where the operator L M is given by 

cos(x) \ _ / cos(x) Dsm(x) 

L l i= \ L c ^ + 2B | Rk,e,n I 

-Dsin(x) / \ cos(x) Dsm(x) 



with 



and 



Here 



L 



Rk,e,^* 

1 „ (X 



l)J,X 



! / X \ 

Rk,eA X ) = R kA X > 6 ) = ^ R [p e ) e 

Thus one has for any eigenvalue A that 

y z(x) J r=0 y z r {x) J 

I w r (x) \ 

The term ( JLl. — A) I is a 2tt periodic function. Since none of the functions share 

\ Z r{x) ) 

a common period shorter than 2irn, the equality 

. . / w r (x) \ 

\ ZrW J 

must hold for each value of r. This shows that one can decompose the spectrum of JL on 
Hl{S 2 nn) as a union of the spectra of the operators JLl posed on i^KS^), i.e. one can 
write 

n-l 

cr(JL) = (J cr(JL r/n ). 

r=0 

Note, one cannot rely on standard Floquet theory since the spectral problem is not an 
ordinary differential equation. In the succeeding sections we study the problem JL^ on S^, 
with ji G [0, 1), in order to deal with arbitrary values of r/n. 



B. Basic Results about the Convolution Kernel and Linearization 

We prove a number of technical lemmas concerning the convolution and linearization that 
are used throughout the remainder of the paper. 
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Lemma 6. Given Hypotheses H2' and HA' , the operator R k ,e,n '■ L^S^) ~~ ► L^S^) 
compact. Further, Rk,e,fi is continuous in \i, \i e [0,1], with respect to the \ |-| | 2 v -norm. 

Proof. Using the same arguments as in Lemma [2j one finds the Fourier series representation 
of Rk,e,nj which we denote as Rk,e,)i, as 

K k,e,n - | 

V 

where A k;€ ,ix is diagonal and (Afc )£j/i )^. = Rk,e(j — /-0- Since in Hypothesis if 2' we assume 
-Rfc j(E G L 1 (R), by the Riemann-Lebesgue lemma [12], we have 



lim (Ajfe,^)^. = 0. 



Defining 



JV 
j=-N 



3\ X )i 



we see that for ||/|| L2(S2w) = 1 



||^ iejM */-< e , M /|| 2 < £ |^e(j-A*)| S 

\lil>iv 

bmce we assume m Hypothesis HA' that |C| < (1 + | s | ) ^, the above sum decays to zero 
as N — > oo, and the operator R k ,e,n* is a uniform limit of finite rank operators. Therefore, 
so is Rk,e,ni an d Rk,e,(i must then be compact. 

To prove the last part of the lemma, we note that for //, fi e [0, 1] 

R k>e (j -y)- R k>e (j - fi) = [°° R k (x; e)e~ ijx [e^ x - e^' x ) dx, 

so using Hypothesis if 2' and the Dominated Convergence Theorem shows that R kte (j — fj,) — > 
Rk,e(j — fJ> ) as /i —>■//, or Rk, e (j — /-0 is continuous in /x. Likewise, we have, using Hypothesis 
HA', 

\\Rk,e,,*-R k ^*\\l<2S(y) + 2S(y) 



where 



^) - X! (1 , ,, , • ^TI7- 



J=-oo 



(1 + \ke(j - y)\Y 
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For ji G [0, 1], one has 



-. oo 

s^) < y ttt- + y — 

VW - (1 + \kei\) 1+2t 1 



so using a dominated convergence argument, one gets | \ Rk,e,^ * — Rk,e,n'* \ | 2 - > as // — »• //. 
Thus | |-Rfc,e lAt — Rk,e,n' 1 1 2 ^ — ^ as — > //, so -Rfc, e ,/i is continuous in /x. □ 

From the previous lemma one gets 

Lemma 7. T7ie operator has a compact resolvent on L\{S2 % )- 

Proof. Given that = L C)M + 2K(e; fi; D), where 

(cos(x) \ _ / cos(x) D sin(aO 
Dsin(x) I \ cos(x) D sin(x) 

one has that K, where we have suppressed the dependence on e and D, is compact since it 
is the product of bounded and compact operators. Further, a straightforward application of 
Fourier series shows that L C)(U has compact resolvent on L 2 (S , 27r ). Let A be a complex number 
with nonzero imaginary part. Then 

/ + {L c ,n - A)" 1 K = (L C , M - A)" 1 - A) . (13) 
The operator / + (L CjA1 — A)~ K is Fredholm since (L CjA1 — A) 1 K is compact. The right- 



hand side of (13) has a trivial kernel since L M is self-adjoint. Thus the left-hand side of (13) 
also has a trivial kernel, which implies that / + {L c ^ — A) 1 K has a bounded inverse [T6] . 
Therefore, from 

(L, - A)" 1 = (/ + (L c ,„ - A)" 1 K) 1 (L c ,„ - A)" 1 , 

one sees that (L M — A) 1 is the product of a bounded and a compact operator, and is therefore 
itself compact. □ 

Assuming that A is in the resolvent of JX M , and using that 

( JL, - A)" 1 = — (Lfj, — 7 r X J (I- in J - A) (L„ - 7 ) _1 J)" 1 , 

where 7 is in the resolvent of L M , one sees that JL M has a compact resolvent on LK^W) 
since (JL M — A) 1 is the product of compact and bounded operators. 

We now need to establish some limiting behavior of the operator Rk,e,n as the nonlocality 
parameter e becomes large. We prove: 
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Lemma 8. Given Hypothesis HA' , for ji ^ 0, Hindoo ||i2fc jejjll *|| 2 = 0. 
Proof. One has 

/oo /»oo 
R k (x;e)e- i V-ri*dx= / C{x)e- ik ^>dx = ((ke(j - 
-oo J —oo 

Examining the L 2 (S , 2 7r ) norm of the operator Rk t e,^*, one gets 

oo 2 

IliWlla ^ E C(* c 0'-^)) • 



j=-oo 



Note, the sum is convergent by Hypothesis HA' . For a given value of the nonlocality param- 
eter e and an arbitrarily chosen value of 5, choose N such that 



E p^'-aO) 

b'l>Jv 



2 S 

<2- 



Next, choose e large enough such that 



N 



j=-N 



2 5 
< -. 

2 



The second assumption does not alter the first since choosing a large e value corresponds to 
choosing a larger value of N. Thus, for (j, ^ 0, 



lim \\Rk,e,n*\\ 2 = 0, 

e— >oo 



and i?fc, e ,/i — >■ uniformly in norm as e — > oo. 



□ 



We finally prove that the resolvents of JL^ and JL converge in the LKSW^norm. This 
is used to show, in effect, that the spectra of one operator is a perturbation in \i of the other. 

Lemma 9. Suppose there exists fi* e (0, 1) such that X is in the resolvent of JL^ for 
< fi < ji* . Further suppose that (JL C)0 — A) -1 exists. Then (JL^ — A)" 1 converges to 
(JL — A)" 1 in the Li^S^-norm as /i — > + . 

Proof. Define the operator -R M (A) = (JL^ — A) -1 . Then we have that 



R^X) - R (X) 
We have that 



< 



2.v 



R (X) 



2,v 



I-(J(L IX -L )R (X) + I 



2,v 



L^ — L — I/ c (/i) + Vi (Rk,t,[i — Rk,efl) V 2 , 
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where L c (fi) = L c ^ — L Cj0 , and V\ and V 2 are such that 2K(e; /i; D) = Vi-Rfc^,^^- Using the 
fact that V\ and V 2 are bounded in L\(S2k) an d Lemma pi 



lim 

/x->0+ 



^1 (Rk,e,n — Rk,e,o) V 2 



0. 



2.v 



Defining i? , c (A) = (JX c ,o — A) 1 , we rewrite ffo(A) so that 



A)(A) = iV(A) / + 2Jtf(e;0;D)i2o, c (A) 



We then get that 

JL c (fi)R (X) 



< 



2.v 



2.v 



I + 2JK(0)R , C (\) 



2,r 



where K(0) = K(e;0;D). The operator L c (/i)_R 0c (A) = L c ([i) (JL c0 — A) -1 is a constant 
coefficient operator. Thus, using the Fourier transform, it is straightforward to show it is 
bounded and must vanish in the ||-|| 2v -norm as fi — > + . Thus we have that 



lim 



Taking /i sufficiently small so that 



J{L^-L )R {X) 



J{L„ - L )R (X) 



2.r 



0. 



< 1, we have that 



2.U 



I-[J(L I1 -L )R (X) + I 



J(L„ - Lq)^o(A) 



< - 

2,v 1 — 



2,1! 



which shows that 



lim 



R^X) - ff (A) 



J(L M - L )R (X) 



0. 



2.v 



2.v 



□ 



IV. STABILITY FOR SMALL POTENTIAL AND LARGE OFFSET SIZE 



A. Computation of the Spectrum with V$ = 



In this section, we compute the spectrum of JL over H 2 (S 27Tn ), with Vq = or D = 1 



(see (12)). As explained earlier, this is done by computing the spectrum of the operators 
JL r / n over iff (-S^), r e {0, • • • , n — 1}. To do this, we notice that we can treat Jf M as a 
constant coefficient operator with a compact perturbation. For the remainder of the section, 
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we assume /i ^ so that the compact perturbation decays uniformly to zero as e — > oo. 
The /i = case is covered by noting that J(L^ — L ) is a relatively compact perturbation of 
JL , which we note was used to prove Lemma [9] Therefore one can find the eigenvalues of 
JLq by taking limits of the eigenvalues of JL^. 

Using the Fourier transform, we compute the spectrum and eigenfunctions of JL C fl ex- 
plicitly. One has 

a(JL c ^) = j±V|(n-/i) 2 -l| :neZ 

and for 0,1, the corresponding eigenfunctions for the eigenvalues on the positive imag- 
inary axis are 

(14) 

while for n — or 1, 

(15) 

Taking conjugates and letting x — > — x gives the corresponding eigenfunctions for the eigen- 
values on the negative imaginary axis. 

The eigenvalue problem for the operator JL^ and eigenvalue X n is of course to find 
nontrivial <p n G -f/f (SW) sucn that 

(JLf, - X n )ip n = 0. 

We write, as in Lemma [7j 

(L Ci/i + 2K(e; 1) + A n J)<p n = 0, 

and let A n = Xoo(n) + X p (n), where Xoo(n) is an eigenvalue of <j(JL CjAt ), and X p {n) is a 
perturbation of X^n) that will be determined exactly. We have 

(L CiM + Xoo(n)J + 2K(q 1) + X p (n) J)<p n = 0. (16) 

Let 

-^^,n = L C)fJ , + A OC) (n) J, 

and 

T(e;n) =2K(e;l) + X p (n)J. 
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Define P n to be the projection onto the null space of L™ n . Since is self adjoint, we 

use a Lyupanov-Schmidt reduction [T7] to rewrite ( Jl6| as 

(17) 
(18) 



where (p 



e„ + M(n)(0 n + ^) = 
P n T(e;n)(0 n + a) = O 

£ n , (p n is in the null space of L™ n , and 



M(n) = (L™ )- 1 (J-P n )T(e;n). 



At this point, the equations (17) and (18) are the same as the original eigenvalue problem. 



No added assumptions or constraints have been made. Therefore solving (17) and (18) is 
equivalent to solving the original eigenvalue problem. 



Rewriting (17) as 



{I + M{n))^ n = -M{n)(f> n , 

we may formally write 

£„ = -(/ + M{n))- x M{n)(j) n = (-M(n) + M 2 (n) )<f> n . 

Though this expansion is valid for sufficiently large e (see Lemma [8]), it is more important 
as a motivation to look at the terms M k (n)(p n . For example, let n = or 1, with A OC) (n) on 



the positive imaginary axis, so that cj) n is given by (15). For D = 1 

cos(x) \ / cos(x) sin(x 

T(e;n)=2B\ ) R k , e ,J 

sin(x) / \ cos(x) sin(x) 




\ p (n)J, 



so that 



where 



T(e; n)(j) n = (Br n+1 - i\ p (n))(j) n + Br n+1 (f) n , 



-i{n+2)x 



and 

r n = Rk,e( n - I 1 )- 

Note, we suppress the parameters e and /z in f n for the sake of clarity in the presentation. 
Thus (/ - P n )T(e; n)<f) n = Bf n+1 <f) n , and 

' k 2 



-((n + 2- nY -l)+zA 0O (n) </> n . 
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Hence, for n = 0, 



and for n — 1, 



M(O)0 O 
M(l)0x 



p(/i-ir ' 
Bf 2 ~ 



P(/i-2) 2 ^ 1 ' 
We consider M(n)4> n for n = or 1. We see that 

T(e; n)<p n = r n+l B<p n + (f n+1 B + i\ p (n))4> n , 

and hence for n = 0, we obtain 

For n = 1, we have 

Define the constants 7„ and S n such that M(n)0 n = 7 n n and M(n)<f) n = 5 n <f) n . Therefore, 
equating 

t ^ n 2 

1 + <)« 



gives a solution of (17). 



Using (18), one obtains 



(T(e;n)0„,0 n ) 



In 



l + 5 r , 



T(e;n) 



0. 



From the work above, since ||0 n || 2 > 0, we see this reduces to 

{Br n+ i - i\ p {n)) — = 0. 



Writing 



Br n+1 + i\ p (n) 

C„(/i) 



where c = k 2 (fi — l) 2 and C\ = k 2 (fi — 2) 2 , we see that we want to solve the quadratic 
equation 

Ap(n) - ic n (fi)X p (n) + c n (fx)Bf n+1 = 0. (19) 
Let X p (n) = X p r \n) + iXp\n), where Xp\n) and X p (n) are real values. Therefore, by 



separating into real and imaginary parts, (19) becomes 

AW(2A«-c n Gu)) =0 



(A«) 2 - (A«) 2 + cMXf + cMBf n+1 = 0. 
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If we assume \ p r \n) ^ 0, we have Xp\n) = — , which implies 

(A^) 2 = -| - cMBf n+1 . 

The right-hand side of the above expression is always negative by construction, since r n+1 > 
0. Thus Xp\n) = 0, so that X p (n) = i\y(n). Therefore we have 

X p( u ) = \ (cn{^) - {c 2 M + Ac n (/i)Bf n+1 ) 1/2 J . 



Since we know, as shown in Lemma |8j that f n — > as e — > oo, and we need \ p (n) — > 
as e — > oo, this determines the correct sign when solving the quadratic equation in Xp\n). 
With this choice of A p (n), we see S n > 0, so 1 + 5 n > 0, and the choice of £ n is well defined. 

For the case that n / or 1, proceeding in a fashion identical to that above, one shows 
that 

M(n)(j) n = 7„0 n , 



and 



Here 



M(n)4> n = 5 n (f> n . 




and 



where c n (/z) = k 2 {n — \i — l) 2 . 



In — f \ i 

_ Br n _ x - i\ p (n) 

C„(/i) 



We again equate = ^(fin, which solves (17), and from (18) one gets a character- 

1 + On 



In 

istic equation for X p (n) which is 

\ 2 p + ic n (n)\ p + c n {n)Br n _ x = 0. 

Finally, 

X p( n ) = \ (-Cn(M) + (Cnifj) +4 Cn (/i) J Bf„_ 1 ) 1/2 j . 

Given that the operator 2K(e; fi; 1) is not symmetric with respect to conjugation followed 
by equating x to —x, we must repeat the above computations except now with the expansions 
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around the eigenvalues along the negative imaginary axis. The process is identical to that 
above, and we only list the results. For n — or 1, we have 



K( n ) = 2 [~ C M + ( c n(ri + 4c n (/i)Sr n _!) J 
where c n (/z) is given by 



' k 2 (fi + l) 2 , n = 0, 



c„(/i) = < 

The corresponding eigenfunctions are 

1 



krfjr, n = 1. 



where 7 n = —^4 and <5 n = 
Likewise, for n ^ 0, 1, we have 




1 / l/2\ 

V n ) = 2 ( c «(^) ~ ( c n(^) + 4c n (ii)Br n+1 ) J 
with c n (n) = k 2 {n — /i + l) 2 . The corresponding eigenfunctions are 




with 7„ = r " +1 and 5 n = . 

The only issue remaining is whether we have captured the entire spectrum of JL^ for 
each value of e e [0, oo). However, every eigenvalue is a perturbation of an eigenvalue in the 
constant coefficient case, which has only simple eigenvalues in its spectrum since JL C)fl is a 
skew-adjoint operator with compact resolvent. Hence we have not missed any eigenvalues 
due to multiplicity. Thus we have computed a^JL^) for V = 0. 



B. Krein Signature 

For a purely imaginary semisimple eigenvalue A G a^JL^) with eigenvector (p, the Krein 
signature of A is defined as sgn((L A1 y9, p)) [18J. Let 

a n {B,[i,e, k) = ln , 
1 + b n 
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which, for the eigenvalues that represent perturbations of eigenvalues on the positive imag- 
inary axis, is given by 

Br n ±i 

u n - 



In I y / 4+4cni?r n±1 

2 ' 2 nil 

where the + in ± corresponds to choosing n — 0, 1 and the — corresponds to n / 0,1. A 
similar expression can be derived for the eigenvalues starting on the negative imaginary axis. 

Given the definition for a n , it is straightforward to show for D = 1 and </? n = <p n + £ n , 
that 

(L^ n , p n ) = 2vr£; 2 ((to - /i) 2 - 1 + a 2 ((n ± 2 - ^) 2 - l)) 

+47r£r n± i(l-a n ) 2 . (20) 

Along the positive imaginary axis, one again lets the — of ± correspond to the case 
to 7^ 0, 1, while we take + for n = 0,1. This relationship is reversed on the negative 
imaginary axis. Thus we see, starting on the positive imaginary axis, for n^0, 1, 2, or 3, 



all the terms in (20) are positive. For n = 2 or 3, we note that < a n < 1, thus 

(to - /i) 2 - 1 + a 2 n ((n - 2 - /i) 2 - l) > 4 (n - fi - 1) , 

which is positive for n = 2 or 3, fi G [0, 1). 

Likewise, along the negative imaginary axis, for n^0, 1, —2, or —1, all terms are positive. 
For n = — 1 or —2, we have 



TO 



u) 2 - 1 + ((to + 2 - u) 2 - 1) > -4 (to - u + 1) 



so that the eigenvalues corresponding to to = —1 and —2 on the negative imaginary axis 
have positive Krein signature. 

However, for to = or 1 on either part of the imaginary axis, if we let B — » + , a n — > 0. 
Therefore we see that for sufficiently small B, with all other parameters fixed, the eigenvalues 
-^oo( n ) + A p (n) for to = or 1 have negative Krein signature. On the other hand, fixing all 
other parameters except B, if we allow the offset size B to become arbitrarily large, then 
a n — )■ 1 and 



lim (to - /i) 2 - 1 + a 2 ((to ± 2 - fi) 2 - l) = 2(to - u ± l) 2 . (21) 
Hence it is possible for eigenvalues to pass through the origin or switch Krein signature. 
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Being more careful, we focus on the eigenvalues with potentially negative Krein signature, 
which are 

A + (l) = i\{2 - /x) (2k - (k 2 (2 - ^) 2 + 4Bf 2 ) 1/2 ) , 



A~(0) = + n) (-2k + (k 2 (l + /i) 2 + 4Bf_i) 
A-(l) = ^(-2fc+(fcV + 4 J Bf ) 1/2 



1/2 



and 



A+( ) = 1 1 <L-A (2 k _ _ + 4B fl ) Va 

One sees that for given k and e, one can find a sufficiently large value of B such that none of 
these four eigenvalues pass through the origin for \i e (0, 1). Noting that r n = Rk,e{n—jj) > 
by Hypothesis H3' , and since Rk,e(n — (jl) is continuous in /i (see Lemma we define 
f n = min^gjo !] Rk, e ( n — /•*)• We define the parameter B* to be 

[3A; 2 3A; 2 k 2 k 2 \ 
B = max < — — , — - — , — , — > . 

^ 4r 2 4r_i r h J 

If B > B*, then all four eigenvalues cannot pass through the origin for /j £ (0, 1). Setting 



/x = |, one has from (21 ) that each of the four eigenvalues must have positive Krein signature 



2> 

for B > B*. 

We now apply a theorem of [19] which states that 

k r + k c + kr = n(L M ), (22) 

where fc r is the number of eigenvalues of JL^ on the positive real axis, k c is the number of 
eigenvalues with real part, k^ is the number of imaginary eigenvalues with negative Krein 
signature, and n{L^ is the number of negative eigenvalues of L^. In order to apply this 
theorem, one needs to show that the operator JX M satisfies Assumptions 2.1a — d in [T9"]. 
Given that JL^ = JL^ C + 2JK, where JK is compact, and that the reciprocals of the 
eigenvalues of JL^ C are square summable, then showing all four assumptions hold for JL^ is 
straightforward. For // = 1/2, k r = k c — k^ = 0, and thus n(Li/ 2 ) = 0. Since the operator 
L M remains invertible for /j e (0, 1), which means no eigenvalue passes through the origin, 
then n(L M ) = for \i G (0, 1). This establishes that every eigenvalue of JL^ has positive 
Krein signature. 
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C. Spectral and Orbital Stability for Small Potential Height 

As shown above, for B sufficiently large, Vq = 0, and /i G (0, 1), there are no eigenvalues 



of negative Krein signature, which by (22) implies that the operator is positive definite. 



Thus a standard perturbation argument guarantees that for small enough Vq, no eigenvalue 



of L M crosses through the origin, and thus we must have n(L M ) = 0. Using (22) again shows 
that ^ is spectrally stable for a given fi with sufficiently small potential height. 

In the case that /i = 0, Vq = 0, one has by continuity of the spectrum with respect to the 
parameter [i that every eigenvalue of JL must be on the imaginary axis. However, for any 
value of Vq, there is an eigenvalue at the origin, with eigenvector 





due to the phase symmetry which generates (j3]). Thus (22) cannot be applied. Likewise, 
there is a generalized eigenvector of JLq at the origin, 

W = 

Using the work above, one formally sees that the eigenvalues at the origin correspond to 
the eigenvalues X{ and \Z\ colliding at the origin for /i = 0. We now prove that the 
generalized kernel of JLq consists only of ip mu and ip gn . First define the projection operator 
(PU, Theorem 6.17) 

p -=2si (ji "- ArldA ' 

where T is a closed, bounded contour in the complex plane such that Tna(JL^) = and the 
origin is inside T. We further suppose that A]~ and XZx are the only eigenvalues of JL^ inside 
T for sufficiently small. Since JLq has a compact resolvent, it has discrete eigenvalues 
that accumulate only at infinity. Therefore, we can also choose V such that V fl a(JL ) = 
and so that T contains only a finite, counting multiplicity, number of the eigenvalues of JL Q . 
Thus the projection Pq is well-defined and finite-dimensional. We then have 



K " P o| | 2 „ < sup || ( JL^ - A)" 1 - ( JL - A)" 1 1 



Aer 



\2,v 



Since 1 1 ( JL^ — A) 1 — ( JL — A) 1 1 1 is continuous in A, on T, which is compact, the supre- 
mum is attained. We further restrict T such that T f] a(JL C:0 ) = 0. Using Lemma [9] then 
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gives 

lim ||P M -P Q || =0. 

/Lt-*0+ ' 

Since P M and Pq are projection operators, we then have ([2T]. pg. 156) that 

dim(Ran(P M )) = dim(Ran(P )), 

where Ran(P At ) denotes the range of P M . By construction dim ( Ran (P M )) = 2, and thus 
dim(Ran(P )) = 2. The dimension of Ran(P ) counts the algebraic multiplicity of an eigen- 
value (see [21], pg- 181), and so we see that the generalized kernel of JL can only consist 
of ip nu and <f gn . 

Since L is self adjoint and has a compact resolvent (see Lemma [7]), it cannot have a 
generalized eigenvalue at the origin for any Vq. At D = 1 (ie. V = 0), it is straightforward 
to show that 

Thus, for Vo = and B > 0, the operator L has a simple eigenvalue at the origin and 
otherwise has only positive eigenvalues. Since the eigenvector at the origin persists for any 
Vo < 0, every nonzero eigenvalue of L remains positive for small values of Vq. This implies 
that n(Lo) = for Vq < and \Vq\ sufficiently small. If n(Lo) = 0, one concludes spectral 
stability for small potential height by way of the following argument. If JLof = \<p, A 7^ 0, 
then 

(L <p,<p) = -A (J<p,<p) , 

and {L ip, ip) > since (p is not in the kernel of L by assumption. (J<p, <p) is strictly 
imaginary and nonzero. Thus A is strictly imaginary. We have now shown that the spectrum 
of JL on D(JX) = iiff (S^Tm) C L^S^Tm), which is decomposed as 

n-l 

a(JL) = (J a(JL r/n ), 

r=0 

is strictly imaginary for small potential height Vo since a{JL r / n ) is strictly imaginary and 
there are a finite number of values r. 

As for orbital stability, again consider (|3]) in the form ip(x,t) = (p ul (x)e~ lult . In other 
words, the solution is generated by the phase symmetry of the Hamiltonian problem ^ 
(see Section 2 for the explicit form of the Hamiltonian). Returning to the original scaling 
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whereby V(x) is a ^ periodic function, with a ^-periodic function, we pose the 



stability problem on Hf(Sx), where T = ^p, with nGN. Thus we are working with (11). 
Then, again, we note that one has the conservation of the quantity 



-r pirn/k 

I(f,9) = v / {f(x)+g 2 (x))dx. 



If we denote the Hamiltonian as %{ipn ipi), one has that W = (p^^ + itj)^ is a critical point 
of E(f, g), where 

E(f,g) = H(f,g)-uI(f,g). 

One has L = "H"(0w,r, — where the primes denote variational derivatives. In the 
given scaling, fi e [0, k), so that a perturbation of period ^p corresponds to = -. In this 
case, we showed that L is positive semi-definite on H$(St) C Lf(Sr). One can forgo the 
requirement that the function d(oj) = E((j) Utr , 4>w,i) be convex (see (TTJ) and conclude orbital 
stability in the space Hi(St) by combining the stable two dimensional real solution into a 
complex function. 



V. SPECTRAL INSTABILITY FOR SMALL OFFSET SIZE 

In contrast to the approach above, we equate the offset size B to zero, for given \i and e. 
We obtain the linearization 

' L+ 



with 




and 



k 2 

L~ = -— ({d x + ifif + l) +2Asin(x)R k ,e^* (sin(x)-). 

We introduced the scaling x — > kx, so that fi e [0, 1). The linearized problem is in "canonical 
form" (see [IS]), and one has the following theorem from 



Theorem 5. fJP} / Let n(L^) and n{L^) denote the number of negative eigenvalues of the 
given operators. With k r the number of real eigenvalues of the canonical system, one has 

h>\n{L+)-n{L-)\. 
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For A small or e large, the problem is spectrally stable since the problem is a small 
perturbation of the constant coefficient case. Further, one has n(L^) = 2. However, with 
e = 0, Lq is a special case of Hill's equation [2U]. We can use the standard spectral theory 
for Hill's equation, from which the spectrum of L~ (i.e. the \i independent operator) is in 
the bands [70, 7j] U [72, 7i] U • • • , where 7, is an eigenvalue of Lq and 7J is an eigenvalue for 
fi = |, with the eigenvalues for all other values of \i filling in the bands continuously If we 
can establish that Lq is positive definite, the same must hold for L~, and we will have then 
shown instability for all values of \x. 

Therefore, setting /x, e = 0, we examine the quadratic form (Lq f, /). Let 

00 



j=-oo 



with 6j(x) from Q where T = 2n. One has 



<w,/>= e ^-^i^r+^i^-^r- (23) 

j=-oo 

There is one negative direction corresponding to the j = mode. Thus, we let / = aeo + 
bg(x), where |a| 2 + \b\ 2 = 1, and g(x) is orthogonal to 6q(x) = It is straightforward to 
show that 

1 f-k 2 _ . . 2 A , A; 2 A 



Therefore, 



L ° 60 = 2^ I ~2~ + 2Asin ^) ) = (— 2" + A)e ° ~ T^ 2 + 6_2) 



A- 2 2/1 
^0 /, f) = M 2 (- T + A) + \b\ 2 (L g, g) - —Re (a*b(g 2 + $_ 3 )) 

^ V Z7T 



and 



A- 2 2/1 
L /, /> > |a| 2 (-- + A) + |6| 2 (L - <?,<?) - -==M|6||«a + $_ 2 |. 

^ V Z7T 



Assume A>^, and define c = — ^ + A We rewrite the above inequality as 

(Lof,f)> (\a\^c--^\b\\g 2 + g^ 2 \\ + \b\ 2 (L g, g) - ^-\b\ 2 \g 2 + g_ 2 \ 



2 



which leads us to examine 



A 2 

L g,g)-—\g 2 + g- 2 \ 2 . 
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With go — 0, we take those terms in (23) that involve only terms in §2 and g~2, which 
reduces our efforts to analyzing 



3k 2 ~\~ A / A 2 1 a i2\ 1 ~ 

1^-2 1 + |^2 1 J - 7^—1^2 + ^-2 



2 



2 Vl " 1 1 7 2ttc' 

Using Young's inequality, \cj2 + <j-2\ 2 < 2 ( | ^ 2 1 2 + \c/2\ 2 )- It follows that if we can satisfy the 

inequality 

3k 2 + A A 2 

—2 — ™ >0 - < 24 > 

we prove that (L^g, g) > 0, and the problem is unstable. A straightforward computation 



shows that (24) is satisfied if 



2 -1 + ,/4-f 

A > -± * = —k 2 « 2.4533A; 2 . 

I — - 

For the sake of presentation we write the instability condition as A > 2. 46k 2 . 

VI. NUMERICS 

In this section, we present numerical results applied to situations for which we expect our 
theoretical results to apply. Then, after calibrating our numerics in this sense, we present 
numerical experiments that correspond to the work in [2]. 



For all the results shown, a filtered pseudo-spectral method [22J is used for the spatial 
variable, while MATLAB's ODE45 function was used for the integration in time. The 
specific filtering function used is a{x) = e ax2j , where a = log(eps), with eps denoting 
machine precision, and 7 = 4. Again see [22] for more details and analysis. 

For the figures in this paper, 128 modes on the domain [0, 8n] were used in the pseudo- 
spectral approximation; higher mode runs were tested and gave identical results to those 
using 128 modes. In each figure, a perturbation of the form 

um(x)e i9ix) 

was added to the initial condition <f>u(x). The function m(x) is a randomly generated, 8tt 
periodic function, normalized so that ||m(x)|| 2 = 1, while v is typically .01. However, in 
certain cases consider v — .1, and where this is the case, it is noted. Finally, given the 
identities derived for the convolution in this paper, the convolution integral turns into a 
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simple term by term multiplication of two vectors in the pseudo-spectral method. Thus no 
approximations to the integral or kernel are made. As in [2J, we convolve against <j>(x) = e~ x2 . 
In every figure, a and k are one. 



Figure la shows the results for B = .01, Vq = —2.46, with nonlocality parameter e = 0. 
As expected from Theorem 4, we see an instability emerge with these parameter values with 
the random perturbation to the initial condition as explained previously. In contrast, Figure 



Tb] shows the case B = 1, Vq = —.01, with nonlocality parameter e = .01, v = .01. The 
numerics behave as the Theorems 2 and 3 predict. We have confidence that the numerical 
results are accurate and correspond to the existing theory. 



Figures [2a| and |2b] show the results for the case B — 1, Vq — —1, e = .01, v = .01 and 
.1. This is a direct comparison to the work in [2j. As can be seen from the figures, the 
underlying solution appears to be robust to perturbations, even with a nonzero nonlocality 
parameter. This contradicts the results of [2J, and seems to imply that ^ is stable in this 
parameter regime. 



VII. CONCLUSION 



We have shown that for a large class of kernels, R(x,e), used to represent long range 
nonlocal interactions in a Gross-Pitaevskii equation, if one lets the range of nonlocality go 
to zero, i.e. e — > 0, then in a rigorous sense, the wavefunction for the nonlocal problem 
approaches the wavefunction for the local problem. This result holds for any smooth peri- 
odic trapping potential V(x). Thus we have demonstrated that generalizing a local model 
to a nonlocal one can be done in a straightforward way, thus expanding the modeling po- 
tential of Gross-Pitaevskii equations. Likewise, we have established the stability properties 
of a particular class of solutions to a nonlocal Gross-Pitaevskii equation. The theory and 
numerical experiments predict that when the offset size B is large these solutions are stable. 
It is therefore possible that under the right conditions these solutions could be observed as 
wavef unctions describing a Bose-Einstein condensate. 
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(a) B = .01, V = -2.46, e = 0, v = .01 (b) B = 1, V = -.01, e = .01, v = .01 
FIG. 1: Confirmation of Theorems 2 and 3 




(a) B = 1, Vo = -1, e = .01, v = .01 (b) B = 1, V = -1, e = .01, u = .1 

FIG. 2: Numerical Predictions for Large Offset Size and Potential Height 
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